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ABSTRACT 

Context. The Galaxy hosts a widespread population of low-energy positrons revealed by successive generations of gamma-ray tele- 
scopes through a bright annihilation emission from the bulge region, with a fainter contribution from the inner disk. The exact origin 
of these particles remains currently unknown. 

Aims. We estimate the contribution to the annihilation signal of positrons generated in the decay of radioactive 26 Al, 56 Ni and 44 Ti. 
Methods. We adapted the GALPROP propagation code to simulate the transport and annihilation of radioactivity positrons in a model 
of our Galaxy. Using plausible source spatial distributions, we explored several possible propagation scenarios to account for the large 
uncertainties on the transport of ~ MeV positrons in the interstellar medium. We then compared the predicted intensity distributions 
to the INTEGRAL/SPI observations. 

Results. We obtain similar intensity distributions with small bulge-to-disk ratios, even for extreme large-scale transport prescriptions. 
At least half of the positrons annihilate close to their sources, even when they are allowed to travel far away. In the high-diffusion, 
ballistic case, up to 40% of them escape the Galaxy. In proportion, this affects bulge positrons more than disk positrons because they 
are injected further off the plane in a tenuous medium, while disk positrons are mostly injected in the dense molecular ring. The 
predicted intensity distributions are fully consistent with the observed longitudinally-extended disk-like emission, but the transport 
scenario cannot be strongly constrained by the current data. 

Conclusions. Nucleosynthesis positrons alone cannot account for the observed annihilation emission in the frame of our model. An 
additional component is needed to explain the strong bulge contribution, and the latter is very likely concentrated in the central regions 
if positrons have initial energies in the 100keV-l MeV range. 

Key words. Astroparticle physics - Gamma rays: ISM - Nuclear reactions, nucleosynthesis, abundances 



1. Introduction 

Over the past decades, it has become clear that our Galaxy can 
produce substantial amounts of antimatter and give rise to sev- 
eral large-scale populations of antiparticles that coexist with our 
matter environment at the current cosmic time. In particular, a 
sizeable number of positrons apparently fill the entire Galaxy 
from the very centre to the peripheral regions. 

The Galactic population of positrons is usually divided 
into low-energy (typically < 10 MeV) and high-energy (typi- 
cally > 100 MeV) particles, and it seems that the re is more than 
histor ical or experimental reasons to do that (Prantzos et al. 
1201 ll) . While high-energy positrons can be observed directly 
as a cosmic-ray component in the interplanetary medium, low- 
energy positrons are revealed indirectly through gamma-ray ob- 
servations of the sky in the ~100keV-lMeV range. The lat- 
ter have indeed shown an unambiguous signature of electron- 
positron annihilation from the inner Galactic regions, in the form 
of a line at 5 1 1 ke V and a continuum below. 

That our Galaxy can produce and host a substantial popu- 
lation of positrons is actually not a surprise. There are many 
physical processes able to provide non-thermal positrons over 
a broad range of energies (photon-photon pair production, ra- 
dioactive decay, hadronic interactions, ...), and as many astro- 
physical phenomena likely to be the sites of one or more of these 
processes (pulsars, X-ray binaries, supernovae, cosmic rays, ...). 
If we add more speculative channels like dark matter particles 
creating positrons through their annihilation or decay, it becomes 



clear that the real challenge of positron astrophysics is not to find 
an origin for these particles, but to identify the source or combi- 
nation or sources that dominates the production. 

In this paper, we focus on the low-energy positron population 
and on one likely source for these particles: the decay of radioac- 
tive species produced by the ongoing nucleosynthesis activity of 
our Galaxy, in particular 26 Al, 56 Ni and 44 Ti. Recent works have 
shown that these elements could account for the annihilation sig- 
nal observed by successive generat ions of gamma-ray telescopes 
jPrantzosl2006tlHigdon et alj2009t) . In this work, we aim at pro- 
viding an additional, alternative assessment of that possibility. 
We adapted the GALPROP public code for cosmic-ray propa- 
gation to simulate the transport and annihilation of radioactivity 
positrons in a model of our Galaxy. Using source spatial pro- 
files based on typical distributions of massive stars and super- 
novae, we explored how the annihilation intensity distributions 
vary upon different prescriptions for the transport. The predicted 
emissions were then compared to the presently available obser- 
vations, coming mostly from the INTEGRAL/SPI instrument. 

We start by recalling the main points of the physics of 
positron annihilation and transport in Sect. and by summa- 
rizing the main observational facts about Galactic low-energy 
positrons in Sect. [3] In Sect. |4] we present the characteristics 
of 26 Al, 56 Ni and 44 Ti in terms of conu'ibution to the Galactic 
population of positrons. We introduce in Sect. [5] the code used 
to model the propagation and annihilation of positrons in the 
Galaxy. Then, we present the simulated cases in Sect.|6]and dis- 
cuss their results in Sect. Q 
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2. Positron physics 

2.1. Positron annihilation 

Galactic low-energy positrons are associated with the strong 
511 ke V celestial signal. For several reasons that will be re- 
viewed in Sect. [3j they are believed to be injected in the Galaxy 
with initial energies below a few MeV. On the other hand, low- 
energy positrons are created non-thermal in all the potential 
source processes, with mean initial energies of a few lOOkeV at 
least. To annihilate, however, most will have to be slowed down 
to near-thermal energies, at least below ~100eV. This energy 
loss proceeds mostly through ionisation/excitation and Coulomb 
interactions in the interstellar gas, a nd it can take ~1 5 yr for a 
typical interstellar density of 1 cirr 3 dJean et al.1 2006). Over that 
time positrons travel at relativistic speeds, which opens the pos- 
sibility for large-scale transport (see below). 

There is a variety of processes by which positrons can an- 
nihilate with electrons: direct annihilation with free or orbital 
electrons, formation of positronium by radiative recombination 
with free electrons or by charge-exchange w ith atoms, ...etc (for 
a complete review, see lGuessoum et al]|2005l) . The physical and 
chemical properties of the medium where the annihilation takes 
place set the dominant annihilation channels and as a result, each 
ISM phase has a typical annihilation spectrum. As we will see 
later, most low-energy positrons seem to annihilate through the 
formation of a positronium atom (Ps), which is a short-lived 
bound state between a positron and an electron. This can oc- 
cur in two different ways: either by charge-exchange, when non- 
thermalised positrons with typical kinetic energy of ~10-100eV 
rip off electrons from H or He atoms; or by radiative recom- 
bination, when thermalised positrons combine with thermalised 
electrons. The Ps atoms occur in 75% of the cases as ortho-Ps 
(parallel particle spins), which decays over /us time scales into 3 
photons producing a continuum below 511 keV, and in 25% of 
the cases as para-Ps (anti-parallel particle spins), which decays 
over ns time scales into 2 photons with 511 keV energies (in the 
centre-of-mass frame). 

Apart from the processes mentioned above, positrons can 
also annihilate directly with electrons at energies higher than 
~100eV and give rise to a Doppler-broadened line extending 
from m e c 2 /2 to Et + m e c 2 /2, where Ek is the positron kinetic en- 
ergy and m e the electron/positron mass. This channel called in- 
flight annihilation is negligible for positron energies of <1 MeV, 
but it becomes increasingly important at higher and higher en- 
ergies. In typical ISM conditions, in-flight annihilation occurs 
for a maximum of ~20-30% of positrons with initial energies 
>10GeV (this prop ortion decreases for stronger magnetic o r ra- 
diation fields, see lSizun et al . 200i|Chemyshov et al1 l2010b . 

2.2. Positron transport 

Once injected in the interstellar medium (ISM), relativistic 
positrons can propagate away from their sources in two ways, 
which we will term ballistic and collisionless transport. In 
the former case, positrons simply follow helicoidal trajecto- 
ries along magnetic field lines and they experience repeated 
interactions with gas particles; in the process, positrons pro- 
gressively lose their energy but experience little deviations^. In 
the latter case, positrons are additionally scattered by magneto- 
hydrodynamic (MHD) perturbations associated with the inter- 

1 See Sect. 3 of I Jean et al.l (2009), the pitch angle of the particles 
remains nearly constant down to ~10keV, and the propagation perpen- 
dicular to the field lines is negligible 



stellar turbulence and they random walk along and across the 
field lines, likely with different properties in the parallel and per- 
pendicular directions; the process can be quite efficient, provided 
there are adequate MHD waves with which positrons can inter- 
act. 

From recent theoretical developments, it seems that inter- 
stellar MHD turbulence can be decomposed into Alfvenic, slow 
and fast magnetosonic modes. These modes are thought to be 
injected at spatial scales of the order of ~ 10 20 - 10 21 cm, as 
a result of differential rotation of the Galactic disk, superbub- 
bles, and stellar winds and explosions. The turbulent energy is 
redistributed to smaller spatial scales through the interaction of 
wave packets in a so-called turbulence cascade. In the process, 
MHD modes can suffer various damping processes that could 
halt the cascade, thus preventing the propagation of MHD modes 
to smaller wavelengths. 

The properties of the turbulence cascade depend on the 
mode, and this has implications for the collisionless trans- 
port of energetic particles. The so-called Kolmogorov scal- 
ing applies to the Alfvenic and slow magn etosonic modes 
dLithwick & Goldreichl 120011: ICho et al.l |2002|) . In that case, 
however, the turbulent energy is preferentially redistributed per- 
pendicularly to the magnetic field, which leads to an ineffi- 
cient scattering of relativistic particles. Conversely, the fast mag- 
netosonic modes f ollow the so-called Kraichnan scaling in an 
isotropic cascade dCho et all 120021: ICho & Lazarianl l2003l). and 
they were shown to have the dom inant contribution to t he scat- 
teri ng of cosmic rays in the ISM dYan & La zarian 2004]). 

iJean et al. (2009, hereafter JGMF09) have made an extensive 
study of both the ballistic and collisionless transport modes for 
positron kinetic energies < 10 MeV, and we briefly summarise 
below their findings about collisionless transport. The authors 
focused on wave-particle resonant interactions occurring when 
the gyroradius of the positron is of order of the parallel wave- 
length of the MHD modes. In fiG fields, this corresponds to 
scales of of ~ 10 9 - 10 10 cm for l-10MeV particles, more than 
10 orders of magnitude smaller than the scales at which turbu- 
lence is injected. 

A minimum wavelength for the MHD modes of the turbu- 
lence cascade is set by Landau damping occurring as the mode 
frequency approaches the cyclotron frequency of thermal pro- 
tons. This corresponds to ~ 10 s - 10 9 cm scales in most ISM 
phases and defines a minimum energy for particle-wave reso- 
nant interactions. In the case of positrons, this minimum energy 
is of order of ~10-100keV, except in hot media with high mag- 
netic field intensities, where it can exceed 1 MeV. Yet, damp- 
ing processes can cut off the turbulent cascade at higher wave- 
lengths if the damping rate exceeds the energy transfer rate, 
and this i s what is thought to happen for some ISM conditions. 
IJGMF09 determined that in the mostly neutral phases of the 
ISM, ion-neutral collisions halt the cascade at spatial scales of 
order~ 10 I6 -10 18 cm for Alfven waves and ~ 10 I7 -10 19 cmfor 
fast magnetosonic waves in warm and cold atomic phases, and 
at scales of order ~ 10 17 - 10 2() cm for both waves in molecular 
phases (where larger scales correspond to smaller gas densities 
for the atomic phases, and to larger cloud sizes for the molecu- 
lar phases). In the ionised phases of the ISM, the Alfven wave 
cascade can proceed undamped down to the minimum spatial 
scales in both the hot and warm media, while the fast modes ex- 
perience significant viscous damping at scales ~ 10 13 - 10 14 cm 
in the warm medium, and strong Landau damping at scales 
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~ 10 17 - 10 18 cm in the hot mediurrQ. 

The above results indicate that the turbulence cascade is 
quenched at spatial scales orders of magnitude greater than those 
required for resonant interactions of 1-10 M eV positro ns, except 
for Alfven waves in ionised medicQ. Yet, JGMF09 argue that 
the scattering by Alfven waves is very likely inefficient due to 
the strong anisotropy of the associated turbulence. The alterna- 
tive of Cerenkov resonance with fast modes, in which no pre- 
cise gyroradius is required for the interaction, could potentially 
take place but under very restrictive conditions involving quasi- 
parallel fast waves and nearly perpendicular positron motion. On 
the other hand, non-resonant interactions with fast modes were 
proved to be q uite efficient at scattering sub-MeV electrons in 
the solar wind and a similar process may well op- 

erate for low-energy positrons in the ISM. Another option for an 
efficient collisionless transport of low-energy positrons is the in- 
jection of turbulence directly at the relevant small spatial scales, 
for instance through the streaming instabilit y but mor e generally 
by any kind of fluid or kinetic instability. JGMF09 concluded 
that positron scattering off cosmic -ray-driven waves is ineffi- 
cient, but that scattering off self -generated waves would deserve 
a detailed investigation. Overall, JGMF09 have presented several 
arguments against the scenario of collisionless transport of low- 
energy positrons, but the issue c annot yet be considered as def- 
initely settled. Independently of JGMF09, Higdon et al. (2009, 
hereafter HLR09) also came to the conclusion that the predomi- 
nantly neutral phases of the ISM do not host the small-scale tur- 
bulence required for the resonant scattering of ~MeV positrons. 
Yet, these authors argued from observations of particle propaga- 
tion in the interplanetary medium that ~MeV positrons should 
diffuse on MHD fluctuations in the ionised phases of the ISM, 
even if the exact nature of the process remains unknown. 

3. Positron observations 

3.1. Constraints on the annihilation sites 

In trying to identify the origin of low-energy positrons, the 
mapping of the annihilation emission by successive genera- 
tions of balloon- and space-borne gamma-ray instruments pro- 
vided a valuable piece of information. The most accurate pic- 
ture available today is provided by the SPI telescope onboard the 
INTEGRAL satellite, which has been in orbit for about 9 years 
and is planned to operate at least until 2014. More than 100 Ms 
of allsky observations concentrated mostly on the Galactic disk 
and bulge are now available. From this, it is clear that the mor- 
phology of the annihilation line emission is dominated by a 
strong central emission consisting of a very peaked part (FWHM 
of 2-3°) at the Galactic centre surrounded by a wider contribu- 
tion (FWHM of 8-10°) from an outer bulge or halo (see Fig. 
[T]l. This inner emission comes on top of a fainter disk-like com- 
ponent extending up to at least |/| < 50°. The most interesting 
feature of this distribution is the relatively high bulge-to-disk 
(B/D) ratio that is obtained for the 511 keV luminosity. Fitting 
projected 3D spatial distributions to the INTEGRAL/SPI data, 

2 Except maybe for quasi-parallel waves, if wave propagation angles 
are not efficiently randomised in wave-wave interactions 

3 Note that Ptus kin et alJ (T2006) investigated turbulence dissipation 
through resonant interactions of cosmic rays, and concluded that this 
process may quench the Kraichnan turbulence cascade of fast modes 
at scales ~ 10 13 cm. While this was invoked as an explanation for the 
~ 1 GeV/nucleon peak in the ratios of secondary-to-primary nuclei ob- 
served in the local cosmic rays, it may also be relevant for the transport 
of lower-energy particles. 



luminosities of 1.2 / 3.1 xlO 43 e + /s for th e bulge, and of 0.8 / 0.5 
Xl0 43 e + /s for the disk were inferred by Weidensnoi ntner et al.l 
(2008b) for two sets of equally-likely modelfl With values of 
order 2-6, the B/D luminosity ratio is higher than those obtained 
for the distributions of classical astrophysical objects or inter- 
stellar gas (with the caveat, however, that the outer disk emission 
is poorly constrained; see Sect. 17.2) . 

Valuable information about where exactly positrons end their 
lives could be obtained from the spectral analysis the 511 ke V 
line, permitted by the high energy resolution of the SPI instru- 
ment. The annihilation emission observed with SPI in the in- 
ner Galactic regions, the 511 keV line and the continuum below 
down to ~400keV, indicates t hat most positrons annihilate in 
the w arm medium of the ISM dChurazov et al.ll2~005l; iJean et ail 
2006). Approximately one half annihilate in the warm ionised 
phase, through the formation of Ps by radiative recombination. 
The other half of the positrons annihilate in the cold and warm 
neutral phases shortly before thermalisation, through the forma- 
tion of Ps by charge-exchange. Most of the positrons that feed 
the 511 keV line emission therefore annihilate through the for- 
mation of a Ps state. The so-called positronium fraction fp s in- 
ferred from SPI is quite high, of the order of >95%, and the re- 
maining few % correspond to direct annihilation with free elec- 
trons in the warm ionised medium. 

3.2. Constraints on the injection energy 

We now turn to constraints on the origin of positrons, and in 
particular on their initial energy and on the reasons why the 
positrons involved in the 511 ke V line emission are generally 
considered to be low-energy particles. For a given positron pop- 
ulation, the level of in-flight radiation at any gamma-ray energy 
> m e c 2 /2 increases with the initial energy of the particles, which 
allows to get constraints on this parameter from observations in 
the 1-100 MeV range. Under the assumption that positrons re- 
main confined to the Galactic bulge reg ions from their injection 
to the ir annihilation, ISizun et all (120061) and Beacom & Yiiksel 
(2006) concluded that the positrons responsible for the 51 1 keV 
emission must be injected with initial energies below a few 
MeV otherwise their in-flight annihilation would give rise to 
an emission excess from the bulge, in contradiction with the 
CGRO/COMPTEL observations in the 1-30 MeV range. A less 
stringent constraint of a few GeV can be obtained assuming an 
injection of positrons ~ 10 5 yrs ago by a non-stationary source 
process, followed by a slowing-down in a magnetic field of a few 
100/l/G in the inner Galactic regions (Chernvsh ov et al.l l2010). 
We will however focus on a stationary source in this work. 

Constraints on the initial energy of positrons annihilating in 
the bulge could also be obtained from other radiations than in- 
flight annihilation. Positrons with energies > 100 MeV would 
emit through inverse-Compton scattering, Bremsstrahlung, and 
synchrotron, and add their contribution to the gamma-ray and 
radio emission from co smic rays interacting with the ISM 
(IStrong et al.ll2007l 1201 lb . Yet, the production rate of positrons 
annihilating at low energy very likely outnumbers that of 
cosmic-ray e lectrons/posi t rons. I n the model of Galactic cosmic 
rays used in IStrong et al.l (120101) . the injection rate of primary 
electrons with energies 100 MeV- 100 GeV is ~ 10 42 e~/s; adding 
secondary electrons/positrons created by inelastic collisions of 
cosmic ray nuclei with interstellar gas approximately doubles 



4 The values are also based on the assumptio n of a positronium frac- 
tion of 0.967, as determined bv lJean et al.l (12006) from the spectral anal- 
ysis of the annihilation emission. 
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Fig. 1. Intensity distribution of the Galactic 511keV emission 
obtained by model-fitting to about 7 years of INTEGRAL/SPI 
observations, using the param eterized components adopted in 
IWeidenspointner et al.l (l2008al) . 



the rate (iPorter et al.l 120081) . As seen previously, however, the 
positron annihilation rate inferred from INTEGRAL/SPI obser- 
vations is a few times 10 43 e + /s; so if these positrons were to 
have initial energies > 100 Me V, they would certainly domi- 
nate the Galactic high-energy or radio emission, at least over 
certain energy ranges. Although current models of the Galactic 
high-energy and radio emission can probably accommodate a 
population of high-energy positrons in addition to conventional 
cosmic rays from supernova remnants, it seems improbable that 
this additional component have similar injection energies and a 
flux an order of magnitude higher. Moreover, the very differ- 
ent morphologies of the 511keV and ~10MeV-100GeV skies 
do not seem to support a high initial energy for positrons, since 
they would first have to lose their high initial energies mostly in 
the disk and ultimately annihilate in a very narrow region in the 
Galactic bulge. 

4. Positron sources 

Positrons are naturally created in the /3 + -decay of some unsta- 
ble proton-rich nuclei. The main species anticipated to signif- 
icantly contribute to the positron production through /3 + -decay 
are 26 Al, 56 Ni and 44 Ti ( 22 Na may also bring_some contri- 
bution, but we will not consider it here; see iPrantzosi 12004). 
Nucleosynthesis is often taken as the most plausible origin for 
the Galactic positrons because the existence and decay of the 
above-mentioned species can be established from various ex- 
perimental results. The most direct evidence comes from ob- 
servations of the characteristic gamma-ray lines or fluorescence 
X-ray lines that accompa ny the decay of the above-mentioned 
isotop es (see for instance Ikenaud et all 120061: ILeising & Shard 
1990). Then, indirect proof can be obtained from the interpre- 
tation of supernovae lightcurves, which are powered by the en- 
ergetic decay products of 56 Ni and 44 Ti, or from the measure- 
ment of present-day cosmic abundances of the daughter nuclei, 
which r esult from the cumulated nucleosynthesis his tory of the 
Galaxy (iMofizuki & Kumagai 120041 : iThe et alJl2006h . The iso- 
tope production yields inferred from these various observations 
would actually be sufficient to provide after f3 + -decay the esti- 
mated ~ 10 43 e + /s that power the diffuse annihilation emission 



observed by INTEGRAL. Yet, as we will see below, these yields 
cannot be straightforwardly translated into a Galactic positron 
injection rate because all positrons do not necessarily escape the 
production sites of their parent radio-isotopes. 

In the following, we present the properties of the positrons 
produced by the decay of the three above-mentioned radionu- 
clides: the plausible large-scale spatial distribution of their injec- 
tion sites, the spectrum they have when they are released into the 
ISM, and their production rate inferred from the estimated nucle- 
osynthesis yields for the parent isotopes. For a mo re detailed dis- 
cussio n about the subject, we refer the reader to iPrantzos et al.l 

<l2oTTh . 



4. 1. The Al radio-isotope 

The 26 Al isotope is believed to be produced predominantly by 
massive stars at various stages of their evolu tion, with typical 
total yields of a few 10" 4 -10" 5 M o per star (IPrantzos & D iehl 
1996). It is released in the ISM by stellar winds and by core- 
collapse supernova explosions. The mean 26 Al lifetime of about 
1 Myr allows it to escape its production site and diffuse ~10- 
lOOpc away before decaying into 26 Mg. The transport of 26 Al 
away from its stellar sources imply that the majority of 26 Al 
positrons are released directly into the ISM and do not suffer 
energy losses prior to their injection. 

The distribution of 26 Al in the Galaxy is truly diffuse and fol- 
lows that of massive stars, as confirmed by the m apping of the 
1809keV decay emission dPluschke et al J 12001) and its corre- 
lation with the microwave free-free emission from HII regions 
or th e infrared emission from heated dust (Knod lseder et al.l 
I1999I) . The Galactic mass distribution of 26 Al inferred from 
the 1809keV emission is strongly concentrated at galactocen- 
tric r adii 3 to 6kpc, at the position of the so-called molecular 
ring dMartin et al.l2009l) . The injection sites of 26 Al positrons are 
thus expected to be distributed in an annular disk and along spi- 
ral arms (such as the thin disk and spiral arms components that 
represent the contribution of localised HII regions surrounding 
massive stars in the NE2001 model for the Galactic distribution 
of free electrons bv lCordes & Lazioll2002l) . 

Because of the long lifetime of 26 Al, positrons are very 
likely injected into the ISM with their original, unaltered /3- 
spectrum. The latter is de scribed by the formula given in 
IChan & Lingenfelterl (1 1 9931 hereafter CL93) and has a mean of 
~0.4 MeV and a maximum of 1 .2 MeV. We note here that most 
26 Al positrons are very likely released in the hot, tenuous, and 
highly turbulent interiors of the superbubbles blown by massive 
star clusters. In these conditions, relativistic positrons may well 
experience stochastic, second-order Fermi acceleration in the 
strong MHD turbulence, or even regular, first-order Fermi ac- 
celeration in the many primary/secondary shocks that travel the 
bubbles. This could significantly modify the initial y8-spectrum 
and bring positrons to higher energies, which could increase 
their range and move their annihilation sites further. This effect 
is, however, beyond the scope the present paper and would de- 
serve a dedicated study. 

The total injection rate of 26 Al positrons can be computed 
from the estim ated total stationary 26 Al mass o f 1.7-2.7 M in 
the Galaxy dMartin et al.ll2009HWang et al.ll2009l) . where the un- 
certainty on the total Galactic 26 Al mass arises from uncertain- 
ties on the exact spatial distribution of the isotope in the Galaxy. 
Taking into account a /3 + -decay branching ratio of 0.82, this 
translates into a 26 Al positron injection rate of (0.20-0.3 l)xl0 43 
e + /s. 
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4.2. The 44 77 radio-isotope 

The 44 Ti isotope is synthesised by explosive Si-burning deep in 
the stellar ejecta during core-collapse and thermonuclear super- 
nova explosions (ccSNe and SNe la respectively), with typical 
yields of a few 1 (T 4 - 1 (T 5 M per event. The mean 44 Ti lifetime is 
about 85 yr, which implies that the radio-isotope remains trapped 
in the young supernova remnant until it decays into 44 Sc and then 
44 Ca. Most 44 Ti positrons are therefore released in the envelope 
of the exploded star and their further transfer to the ISM de- 
pends on the unknown transport conditions in the ejecta and on 
the explosion properties. This must be taken into account when 
translating 44 Ti yields into positron injection rates (see below). 

The injection sites of 44 Ti positrons sample the time- 
averaged distribution of supernova explosions in our Galaxy. 
Core-collapse supernovae proceed from massive stars and are 
thus distributed like them, in a thin annular disk and along spiral 
arms (see 14. II . Thermonuclear supernovae have a slightly more 
subtle evolutionary pathway and actually depend on both the 
old and young stellar populations. Sull ivan et al.l (120061) showed 
from observations of external galaxies that the rate of SNe la de- 
pends on both the star formation rate and the total stellar mass^]. 
The spatial distribution of SNe la is therefore a combination of 
that of massive stars and that of stellar mass. In our Galaxy, most 
of the stellar mass is distributed in an exponential disk with a 
central hole and in a ellipsoidal bulge, with approximately the 
same mass in each component: 2.15 and 2.03 x 10 I() M , respec- 
tively (see the determination of t he three-dimensional shapes and 
parameters in Robi n et al.l2003l) . Overall, the source distribution 
for 44 Ti positrons has three components: a thin annular disk (for 
ccSNe and prompt SNe la), and an exponential disk with cen- 
tral hole and a bulge (for delayed SNe la). Further down, we 
quantify the relative contribution of each component to the total 
source term. 

Because of the intermediate lifetime of 44 Ti, the population 
of decay positrons entering t he ISM may be affe cted by the es- 
cape from the stellar ejecta. Ma rtin et al.l (1201 Ol) computed the 
escape or survival fractions of 44 Ti positrons for two extreme 
transport modes in the ejecta: either positrons can freely stream 
through the stellar envelope, or they are trapped at their birth 
place by some strong magnetic turbulence. In the free stream- 
ing case, the escape fractions are quite high, ranging from 97% 
for light 2M Q ejecta to 83% for more massive 14M ejecta; 
in the trapped case, the survival fractions are lower, ranging 
from 91% to 36% for the same mass range (assuming in both 
cases a typical 10 51 erg explosion kinetic energy Although 
44 Ti positrons do suffer energy losses on their way out of the 
stellar ejecta, the mean energy of the escaping/surviving pop- 
ulation remains close to the initial one, whatever the transport 



5 A scenario based only on stellar mass, that is on the older and lower- 
mass star population, is ruled out at >99% confidence level. 

6 These estima tes for the escape fractions are consistent with the re- 
sults obtained by ICL931 except for the 14 M Q ejecta mass, for which 
they found a survival fraction of 8% only in the trapped case. The ori- 
gin of the discrep ancy may lie in the approximation they used in their 
calculation, while M artin et"aT] (|2010) performed the complete integra- 
tion over time and energy. Another difference between both studies is 
that|CL93 added a so-called slow positron survival fraction for positrons 
that are thermalised but do not annihilate in the continuously-thinning 
ejecta. Yet, in the ejecta assumed to be neutral, the dominant annihi- 
lation process is positronium formation, which occurs over compara- 
tively very short time scales once positrons have been slowed down 
below ~100eV. So if neutral atoms are indeed the dominant species 
in the rapidly-cooling ejecta, no more than a few % of the thermalised 
positrons are expected to survive. 



mode. So most 44 Ti positrons are injected into the ISM with a 
spectrum close to the original, unaltered /^-spectrum, which has 
a mean of ~0.6MeV and a maximum of 1.5 MeV. We note, how- 
ever, that [Zirakashv ili & Aharonianl (1201 lb considered the pos- 
sibility that positrons from the 44 Ti decay chain constitute a pool 
of mildly relativistic particles that can be accelerated to ultra- 
relativistic energies in the remnants of supernova explosions of 
all types. In their model, the ~1 MeV positrons from 44 Sc de- 
cay undergo stochastic pre-acceleration up to ~ 100 MeV in the 
turbulent upstream region of the reverse shock, and are then fur- 
ther accelerated up to multi-TeV energies by the diffusive shock 
acceleration mechanism at the shock. This effect is, however, be- 
yond the scope the present paper. 

The total injection rate of 44 Ti positrons can be com- 
puted from the frequencies of ccSNe and SNe la, and the 44 Ti 
yields and positron escape fractions for each type of event. 
From Tammann et al. (1994), ccSNe are estimated to occur at 
a rate ~2.1 ccSNe/century while SNe la occur at a rate ~0.4 
SNe l a/century. Using the empirical relation from lSullivan et all 
on which we based the spatial distribution of SNe la, the 
occurrence rate of SNe la is 



#SNia = M* x (5.3 ± l.l).10- I4 yr- 1 M Q 1 

+SFR x (3.9 ± O.7).lO _4 yr _1 (M yr 1 )" 1 (1) 



we also find a rate of ~0.4 SNe la/c entury if we tak e a total 
star formation rate S FR = 4M Q yr" 1 dPiehl et al.ll2006l) and the 
above-mentioned to tal stellar mass M* = (2.15+2.03)x 10 10 M Q 
dRobin et alJ 120031) . Then, a star formation rate of 4M yr~ 1 
corresponds to ~2.1 cc SNe/century for a typical initial mass 
function (see Table 1 in iDiehl et alJl20 06). Regarding the 44 Ti 
yields of supernovae, observational estimates are very scarce 
(two direct detectio ns, Cassiopeia A and Gl .9+0.3, and one indi- 
rect, SN1987A: see iRenaud et alJl2006t iBorkowski et al.ll2Q10b 
iMotizuki & Kumagail 120041) . so it is risky to compute some 
average Galactic 44 Ti production rate by integrating theoreti- 
cal yields over a range of supernova progenitors. Instead, we 
simply assumed that ccSNe and SNe la eject on average 2.0 
xlO~ 4 M and 2.0 xlO~ 5 M per event, respectively. These val- 
ue s agree with t he compilation of observational constraints done 
in iMartin et al.l (1201 0[) and with the recently-estimated range 
for Gl .9+0.3 feorkows ki et al.ll2010b . Overall, this corresponds 
to a mean Galactic 44 Ti production rate of 4.2 xlO" 6 M yr" 1 , 
roughly consis tent with the estimate of 5.5 xlO" 6 M yr _1 by 
iThe et all (2006) based on the present-day solar 44 Ca abundance 
and a Galactic chemical evolution model. This 44 Ti production 
rate can be translated into a positron production rate by apply- 
ing a factor 0.94 for the /3 + -decay branching ratio of the decay 
chain and using an escape fraction of 100%. In the absence of 
strong observational constraints, the latter value was chosen be- 
cause the escape fraction in models is >50% for most ejecta 
masses, explosion energies, and transport conditions, so using 
100% gives the correct order of magnitude and allows an easy 
scaling of the results. Eventually, the 44 Ti positron injection rate 
is 0.34 x 10 43 e + /s (with only -2% being contributed by SNe 
la), and we consider an uncertainty range of ±50% due to the 
uncertainty on the Galactic 44 Ti production rateQ. 



7 In The et al. ( 2006), the authors found a ~ 2 times larger uncertainty 
range for the Galactic 44 Ti production rate due to unknown parameters 
in the Galactic chemical evolution model; they also note that the en- 
tire range may be shifted to lower values if there is a nucleosynthesis 
channel to produce 44 C a directly. 
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Table 1. Bulge, disk, and total annihilation radiation luminosities for 56 Ni, Ti, and ' 6 A1 positrons in the three transport con- 
figurations tested. Also indicated are the corresponding bulge-to-disk ratios, and the fractions of positrons that annihilate in the 
Galaxy. 



Source 


transport 


Bulge (K < 3 kpc) 


Disk (K > 3 kpc) 


total ualaxy 


Bulge/Disk ratio 


Annihilation traction 




A 


0.32 


0.58 


0.91 


0.56 


1.00 


56 Ni only 


B 


0.16 


0.47 


0.63 


0.35 


0.68 




C 


0.09 


0.39 


0.48 


0.22 


0.57 




A 


0.01 


0.18 


0.20 


0.08 


1.00 


^Ti only 


B 


0.01 


0.11 


0.12 


0.10 


0.63 




/ ■ 

L. 


U.U1 


U.iU 


U.l i 


U.Uo 


U.jy 




A 


0.01 


0.13 


0.14 


0.08 


1.00 


26 Al only 


B 


0.01 


0.10 


0.11 


0.09 


0.77 




C 


0.01 


0.09 


0.10 


0.07 


0.73 




A 


0.35 


0.90 


1.25 


0.39 


1.00 


56 Ni+ 44 Ti+ 26 A , 


B 


0.18 


0.68 


0.86 


0.27 


0.69 




C 


0.10 


0.59 


0.69 


0.17 


0.56 



Note to the table: The luminosities in columns 3 to 5 are given in 10 43 ph s 1 for the mean positron injection rates determined in Sect. |4] They 
correspond to 51 1 keV emission from parapositronium and direct annihilation, for a positronium fraction //>,=0.95. 



4.3. The 56 Ni radio-isotope 

The 56 Ni isotope is synthesised by explosive Si-burning deep in 
the stellar ejecta during core-collapse and thermonuclear super- 
nova explosions, with typical yields of the order of 10~' M per 
event. The characteristic time of the decay chain to 56 Co and then 
to 56 Fe is <1 yr, which implies that all positrons are released in 
the stellar ejecta at the late supernova / early remnant phase. This 
makes escape a critical point, which actually counterbalances the 
large 56 Ni yields. 

Although 56 Ni is produced by both ccSNe and SNe la, we 
will see below that SNe la are expected to dominate over ccSNe 
in terms of contribution to the Galactic positron population 
through that isotope. Consequently, the injection sites of 56 Ni 
positrons follow the time-averaged distribution of SNe la in our 
Galaxy. The source distribution has three components like that 
of 44 Ti positrons - a thin annular disk, an exponential disk, and 
a bulge - but the relative contribution of each component to the 
total positron production rate is different, as we will see below. 

Because of the short timescale of the 56 Ni decay chain, the 
population of decay positrons eventually entering the ISM is 
strongly affected by the travel through the stellar ejecta. The es- 
cape of 56 Ni positrons is favored by ejecta mixing, which lifts 
iron up to the lower-density ejecta surface, and by low con- 
finement, which decreases the column density experienced by 
positrons in the ejecta. Models exploring these effects result in 
escape fractions up to ~10% for SNe la, while escape fractions 
for ccSNe reach a few % only for the light and rare type Ib/Ic 
explosions a nd in t he limiting case of a fully mixed ejecta (see 
for instance ICL93h . The comparison of such predictions with 
the late lightcurv es of SNe la ind icates an escape fraction in 
the range ~2-6% dMilne et al . 1999, bu t see also the warning of 
lLairetal.1 (120061) ). Alternatively. iMartin etail (1201 0b translated 
the non-detection of 5 1 1 keV emission from the youngest and 
most nearby supernova remnants into an upper limit on the es- 
cape fraction of 12% for SNe la. Overall, SNe la are thought 
to be the dominant source of 56 Ni positrons when compared to 



ccSNe: their lower occurrence rate (0.4 versus 2.1 SNe/century) 
is compensated by their higher iron yield per event (0.6 versus 
0. 1 M , typically), and their average positron escape fraction is 
very likely one order of magnitude higher. Due to the strong en- 
ergy losses experienced in the fresh and dense ejecta of SNe la, 
escaping/surviving 56 Ni positrons very likely have a spectral dis- 
tribution that differs from the original yS-spectrum. Depending 
on the actual mixing and density profile of the ejecta, the mean 
kinetic energy of the positron population may be shifted from 
~0.6MeV at decay to ~0. 3MeV or even below when entering 
the ISM (see the Fig. 4 of lCL93h . We will discuss later on the 
impact of this effect on the predicted annihilation emission. 

With the previously-adopted SNIa occurrence rate of 0.4 per 
century and a typical 56 Ni yields of 0.6 M , the total positron 
injection rate can be computed by taking into account a factor 
of 0.19 for the y6 + -decay branching ratio of the decay chain and 
assuming an escape fraction of 5%. This corresponds to a 56 Ni 
positron injection rate is 1 .53 x 10 43 e + /s. Yet, the positron escape 
fraction is a crucial parameter that is not strongly constrained by 
observations. We therefore considered an uncertainty range of 
one order of magnitude for this factor, for a range of 1 to 10%. 
This translates into an uncertainty range on the 56 Ni positron in- 
jection rate of (0.31 -3.10) x 10 43 e + /s. 

5. The transport code 

The propagation of positrons in the Galaxy was modelled with a 
modified version of the publicly-available GALPROP code (see 
http://galprop.stanford.edu). The code was originally created to 
analyse in a consistent way the growing body of data available 
on Galactic cosmic rays (direct particle measurements, diffuse 
emissions in radio and at high energies). Nevertheless, it is gen- 
eral enough to be adapted to other studies of particle transport in 
the Galaxy. 

In the following, we briefly review the constituents of the 
GALPROP model that are relevant to the study of nucleosynthe- 
sis positrons. We also introduce the modifications implemented 
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to simulate the propagation and annihilation of positrons in the 
Galaxy. 

5. 1 . Transport processes 

GALPROP numerically solves a general diffusion-convection 
equation for a given source distribution. The transport mecha- 
nisms included in the code are diffusion in position and momen- 
tum, convection away from the plane and the associated adia- 
batic cooling, momentum losses from many processes, radioac- 
tive decay, and nuclear fragmentation. The equation is solved us- 
ing a Crank-Nicholson implicit second-order scheme, with free 
particle escape assumed at the spatial boundaries. A typical sim- 
ulation is run for decreasing time steps until a steady state is 
achieved for all species over the entire computational grid. The 
present modelling of positron propagation in the Galaxy includes 
diffusion in position and energy losses, but no convection or dif- 
fusive reacceleration (and of course no nuclear process). 

By construction, the transport in position space is treated as 
a diffusive process and we will explain in Sect. !6.2l how this can 
be justified for the cases we considered. We modified the code 
to allow the simulation of inhomogeneous diffusion (with dif- 
ferent properties in the bulge and in the disk for instance). To 
avoid potential problems arising from discontinuities in the grid 
of diffusion coefficients, we implemented smooth transitions be- 
tween regions with different diffusion properties. In the case of 
a more efficient positron diffusion in the disk compared to the 
bulge (see Sect. I6.2I ). the diffusion coefficient is described by a 
Gaussian-type function having its lowest value at the Galactic 
centre and reaching its highest value beyond a certain scale ra- 
dius and height above the plane. 

Regarding the transport in momentum space, GALPROP 
includes the main energy loss processes for high-energy 
charged particles, such as Bremsstrahlung, inverse-Compton, 
and synchrotron, but the dominant one in the keV-MeV range 
are Coulomb interactions and ionisation/excitation of atoms. 
Io nisation/exc i tation are implemented following the prescription 
of Pa ges et all (Il972h . based on Bethe's theoretical formula, with 
experimental values for the ionisation potentials of neutral H and 
He, and without correction for the density effec t. Coulomb losses 
are implemented following the prescription of Ginzb urgl (119791) 
for the cold plasma limit. As noted below, we also implemented 
direct in-flight annihilation in the code. This constitutes an ad- 
ditional catastrophic loss process, but its contribution for MeV 
positrons is insignificant. 

Overall, the transport equation solved for positron propag a- 
tion in the Galaxy (until a steady state is achieved) is the follow- 
ing 



(2) 



where <p is the positron distribution function, D is the spatial dif- 
fusion coefficient, E p is the positron energy loss rate, and Q is 
the source term. All these quantities depend on position r and 
positron energy E p . The equation does not contain an explicit 
term for annihilation. As we will see below in Sect. 15.31 the 
annihilation rate is computed at each position from the rate of 
positrons being slowed down below ~100eV. 

5.2. Interstellar medium 

GALPROP includes average analytical spatial distributions for 
the main gas states: molecular (H2), atomic (HI), and ionised 
(HII). These distributions are used in the computation of the 



propagation of cosmic -rays throughout the Galaxy, for the de- 
termination of energy losses for instance. Then, for the predic- 
tion of diffuse emissions from cosmic rays interacting with in- 
terstellar gas (such as Bremsstrahlung or n° production and de- 
cay), GALPROP allows to recover the fine spatial structure of 
the gas through the use of the observed gas column densities (in 
the HI 21 cm and CO 2.6 mm emission lines). In our simulations 
of the annihilation emission, however, we used only the analyti- 
cal gas distributions because the observational constraints at our 
disposal, coming mostly from INTEGRAL/SPI, have an angular 
resolution several times above that of HI or CO surveys. 

Atomic hydrogen is the dominant gas phase in terms of to- 
tal mass and has a relatively large filling factor. The GALPROP 
code uses a 2D ana lytical distribution for H I. The radial distri- 
bution is taken fro mlGordon & Burtonld 19761). while the vertical 
distribution is from Dick ey & Lockmanl(ll990l) for < R < 8 kpc 
and ICox et all (1 19861) for R > 10 kpc. with linear interpolation 
between the two ranges. 

Molecular gas accounts for about a third of the total gas 
mass of the Milky Way, but it is concentrated in complexes of 
dense, massive clouds with low filling factors. The GALPROP 
code emplo ys a 2D analy t ical d istribution for H2, using the 
mode l from iFerriere et al. (120071) for R < 1.5 kpc, the model 
from lBronfman et al.1 (1 1988b for 1.5 kpc < R < 10 kpc, and the 
model from lWouterloot et all (1 19901) for R > 10 kpc. 

The ionised hydrogen makes up about 10% of the total 
gas mass of the Milky Way and is therefore the least massive 
component. Yet, It actually occupies most of the Galactic vol- 
ume and has a very extended distribution across the plane. The 
GALPROP code uses a two-component model for the HII dis- 
tribution: a thin disk tracing localised HII regions in the plane, 
mostly in the molecular ring and spiral arms, and a thick disk 
representing the more diffuse warm ionised gas that exists out- 
side the well-defined HII regions. The former is modelled by 
the 2D/axisymetric thi n disk component of the NE2001 model 
(ICordes & Lazio 2002), while the latter is modelled by a verti- 
cal exponentia l dist ribution with the scale height determined by 
iGaensleret ai1(l2008l) . 

5.3. Annihilation 

In its publicly-available version, GALPROP does not include an- 
nihilation processes. As explained in Sect. 12. II positrons can an- 
nihilate with electrons through a variety of processes, the rel- 
ative importance of which depend on the characteristics of the 
medium where the annihilation takes place. Yet, GALPROP 
does not include a fine description of the ISM at the typical 
scales of the various gas phases, but only the averaged and ax- 
isymetric gas distributions described above. At a given position 
in the space grid, the annihilation of positrons could not occur 
in a definite, well-identified phase but in an average medium re- 
sulting from the superposition of the large-scale distributions of 
molecular, atomic and ionised gas. 

We implemented positron annihilation in the following way: 
the rate of positron annihilation at a given position in the Galaxy 
is based on the rate of positrons being slowed down below 
~100eV, that is the energy below which positrons may expe- 
rience charge exchange in-flight. Theoretically, this applies only 
to the neutral phases because in ionised phases positrons form 
Ps by radiative recombination or annihilate directly after com- 
plete thermalisation. In the warm ionised phase, the thermalisa- 
tion from 100 eV energies and the subsequent annihilation occur 
on short time scales compared to the slowing-down, so limiting 
the modelling to >100eV energies has no consequences on our 
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results. In the hot phase, however, the thermal or near- thermal 
positrons have long lifetimes owing to the very low density. In 
the disk, these positrons are thus expected to be transported out 
of the hot cavities, very likely by strong turbulence, and to anni - 
hilate in the denser surrounding phases dJean et al ] |2006ll2009h . 
This happens on time scales < 1 Myr, of the order of the slow- 
ing down time, and so once again limiting ourselves to >100eV 
energies will not strongly alter the results (especially since the 
cell size in our spatial grid is larger than the typical size of hot 
cavities, ~100pc). 

The corresponding positron annihilation emissivity q for po- 
sition r and photon energy E y reads 

9(r,E 7 ) = [v(r,E p )E p (r,E p )] mey (3) 

where fp s is the Ps fraction, and the s(E y , fp s ) function gives the 
spectral distribution of the radiation. Yet, since we do not track 
the annihilation in individual ISM phases, we cannot directly 
model the annihilation spectrum over the Galaxy. We therefore 
fixed the spectral characteristics of the annihilation radiation: the 
Ps fraction was taken to be //>,=0.95, the 511keV line was as- 
sumed to have a Gaussian shape a nd a width of 6 ke V, and the Ps 
continuum was described by the lOre & Powelll (1 19491) formula. 
We then focused on the intensity distribution over the sky, which 
is obtained by integrating the above emissivity along the line-of- 
sight for each direction to the sk>0. 

We also implemented in our version of the GALPROP code 
the process of direct in-flight annihilation. In that case, it was im- 
plemented with its complete differential cross-section because it 
depends only on the bound and free electron density. However, 
we do not insist on that aspect because in-flight annihilation is 
negligible for Me V positrons. 

6. Simulation setups 

The numerical model described above allowed us to simulate the 
propagation of nucleosynthesis positrons and compute their an- 
nihilation emission for the entire Galactic system in a clear and 
consistent way. The predicted annihilation emissions basically 
rely on two inputs: the source spatial distribution and the prop- 
agation parameters. The source distributions were considered to 
be more firmly established than the transport conditions, so we 
fixed the source spatial parameters to the values given in Sect. [4] 
and 16.11 and explored how the annihilation emission vary upon 
different prescriptions for the transport. 

All simulations presented in this paper were made in 2D 
cylindrical geometry. The spatial grid extends from to 20 kpc 
in Galactocentric radius r, with a step size of 250 pc, and from -4 
to 4kpc in Galactic height z, with a step size of 100 pc. This im- 
plies that the Galactic halo extends up to 4 kpc on either side on 
the Galactic plane. The spatial boundary conditions at the edges 
of the spatial domain assume free particle escape, which means 
that the ultimate fate of positrons is either to annihilate in the 
disk or halo if they are slowed down below lOOeV, or to escape 
if they diffuse far enough. The resolution of our space grid cor- 
responds to ~l-2° at the distance of the Galactic centre, which 
is about the size of the strongest peak in the 511 ke V signal. The 
energy grid runs from 100 eV to 2MeV on a logarithmic scale, 
with 156 energy bins. 

8 Note that our results can be rescaled to any other Ps fraction than 
the /p^O.95 adopted: for instance, the 511 keV intensities presented 
here can be converted to the case of a different Ps fraction g Ps after a 
multiplication by (l-0.75g Pj )/(l-0.75/p s ). 



6.1. Sources parameters 

In this section, we present the 2D source spatial distributions that 
were adopted for the different positron source. We saw in Sect. 
[4] that the injection sites of radioactivity positrons follow a spa- 
tial distribution in the Galaxy that has three components: a star- 
forming disk consisting of an annulus and spiral arms (for parent 
isotopes coming from massive stars, their ccSNe, and also from 
prompt SNe la), together with an exponential disk with central 
hole and an ellipsoidal bulge containing most of the stellar mass 
(for parent isotopes coming from delayed SNe la). 

For the star-forming disk, we used the azimuthally-averaged 
radial profile of the Galact ic star formation rate determined by 
Bois sier & Prantzosl (I1999I) : since massive stars have short life- 
times of a few Myr, their distribution and that of ccSNe closely 
follows that of the star forming activity. The vertical profile was 
assumed to be exponential with a scale height of 200 pc. 

The stellar mass disk and bulge components were modelled 
by the functions: 

f D (r, z) = «£>,() X (e~*°* - j x e' 7 ^ (4) 

fair < R B ,c, z) = «b,o X e ^* z "'< 1 (5) 
f B (r > R B ,c, z) = n By0 x/'Wx e -W- R »,c) 2 (6 ) 

where (r, z) are the Galactocentric distance and the height above 
the Galactic plane and the indices D/B stand for disk and bulge, 
respectively. These formula approximate the usual prescriptions 
for th e stellar disk and bulge shapes (see for instance Robi n et alj 
2003). Each component D/B is described by 1 density normal- 
isation factor no and 3 geometrical parameters (Ri,,Zi,,R c ): a 
scale radius and height, a cutoff radius. We adopted the follow- 
ing triplets: (2.5,0.3,1.3) for the disk, and (1.6,0.4,2.5) for the 
bulge, in units of kpc. The normalisation factors are computed 
from the data given in Sect. |4]for each parent radio-isotope. 

6.2. Diffusion parameters 

The recent theoretical studies reviewed in Sect. 12.21 could not 
conclusively determine if the transport of low-energy positrons 
is ballistic or collisionless. Strong arguments were presented in 
favour of a ballistic transport, in particular in the neutral phases, 
but these conclusions are not yet backed by solid experimental 
evidence. We therefore tested two extreme transport configura- 
tions in the diffusion approximation, with small and large coef- 
ficients corresponding to collisionless and ballistic transport, re- 
spectively. This allowed us to assess the impact on the resulting 
annihilation emission. In addition, we simulated an intermediate 
case of inhomogeneous diffusion, where the transport was as- 
sumed to be collisionless in the Galactic bulge and ballistic in 
the Galactic disk. Below, we present these transport scenarios in 
more detail, including the diffusion coefficients used in the code. 
We recall here that GALPROP does not include a fine descrip- 
tion of the ISM with distinct phases, and so diffusion should be 
implemented with average properties over a few lOOpc scales. 
In addition, we assumed in all cases an isotropic diffusion, but 
anisotropies may well exist for diffusion along the Galactic plane 
and off that plane towards the halo, or along spiral arms and 
across them. 

Case A: We tested the scenario of a collisionless transport 
with homogeneous properties over the whole Galaxy. For the 
ionised phases of the ISM, this option actually remains partly 
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open in the work o f|jGMF09 and is invoked in lHLR09l (see Sect. 
I2.21 >. As a limiting case, we assumed that it holds for all ISM 
phases over the entire Galaxy. In this scenario, we implicitly as- 
sumed the existence everywhere in the Galaxy of MHD turbu- 
lence with a sufficient energy density in the required wavelength 
range, so that low-energy positrons can be efficiently scattered. 
The corresponding random walk process was then treated in the 
diffusion approximation with a coefficient of the form: 

D w (R)=/3D (^-\ (7) 

D W (R) = 3.7 x KFcmV 1 ) (8) 

where R = p/e is the rigidity of the particle and f3 = v/c. The 
adopted normalisation factor Do as well as the spectral index 6 
are taken from cosmic-ray propagation studies interpreting lo- 
cal cosmic -ray measurements at ~l-100GeV in the frame of 
a Galactic cosmic-ray propagation model. The spectral depen- 
dence of the diffusion coefficient is assumed to hold down to 
~MeV energies, which so far remains unproven (see Sect. 12.2) . 
On the other hand, the resulting diffusion coefficient is small 
enough that the typical range of a MeV positron in a 1 cm -3 
density medium is of order ~10pc, meaning that nucleosynthe- 
sis positrons annihilate close to their sources (technically, they 
hardly escape the cell where they were injected). We assumed 
that diffusion is homogeneous over the Galactic volume, which 
is equivalent to assuming that the characteristics of magnetic tur- 
bulence are homogeneous over the Galaxy. 

Case B: We then tested an intermediate, inhomogeneous 
transport scenario, where the transport is collisionless in the 
Galactic bulge and ballistic elsewhere, notably in the Galactic 
disk. The limit between the two regions is set at a Galactocentric 
radius of about 3 kpc. The form of the diffusion coefficient in the 
ballistic case is given below. This scenario actually is motivated 
by the observation that the bulge and the disk have different ISM 
compositions: the bulge is largely dominated by the hot ionised 
phase, while the disk is mostly filled with neutral gas. Estimates 
for the filling factors of the warm neutral medium, warm ionised 
medium, and hot medium are respectively (0.2,0.1,0.7) for the 
bulge and (0.5,0.3,0.2) for the disk, while the molecular medium 
and c old neutral medium occupy compara tively negligible vol- 
umes (iJean et alJl2006l:lHigdon et alJl2009h . From Sect.O this 
difference may imply that positrons would be more efficiently 
confined in the bulge by small-scale MHD waves in the predom- 
inant ionised phases, while they would stream more easily in the 
disk where damping mechanisms prevent small-scale MHD per- 
turbations in the predominant neutral phases. In reality, large- 
scale diffusion properties over a few lOOpc scales reflect the 
actual ISM composition in each Galactic region (like the fact 
that the ionised phases of the disk occupy 50% of the volume 
and thus would moderate the easier streaming experienced in 
the neutral phases). In our simulations, however, we pushed the 
contrast between bulge and disk to the maximum and set the for- 
mer to be fully collisionless and the latter to be fully ballistic. 
This choice clearly is biased towards reproducing the high ob- 
served bulge-to-disk 511 keV luminosity ratio, but we will see 
that even with that assumption, our predictions do not agree with 
the INTEGRAL/SPI results. 

Case C: Last, we tested the scenario of a ballistic transport 
with homogenous propert ies over th e whole Galaxy. This case is 
suggested by the work of JGMF09, who concluded that ~MeV 
positrons are likely not scattered by MHD turbulence whatever 
the ISM phase. In this situation, ~MeV positrons have large 



mean free paths of several kpc in terms of pitch angle scattering. 
On small scales, the transport is anisotropic and positrons follow 
magnetic field lines in a ballistic motion. On larger scales, how- 
ever, we expect an isotropization of the transport by strong ran- 
dom fluctuations of the Galactic magnetic field on typical scales 
~100pc (an argument t hat is commonly in voked in cosmic -ray 
studies; see for instance Ptus kin et al.l l2006). In the ballistic sce- 
nario, we therefore approximate the transport by a diffusion with 
a characteristic mean free path Lb, the largest scale of magnetic 
fluctuations. The coefficient then takes the form 

D C (R) = \pch B (9) 

D C (R) = 3.1 x 10 30 cm 2 s _I x/?|-^-| (10) 

\100pc / 

In this case, the diffusion coefficient has a very modest rigid- 
ity/energy dependence down to a few lOkeV, and the range of 
particles is therefore imposed by the energy dependence of the 
slowing-down processes. We assumed that diffusion is homoge- 
neous over the Galactic volume, which is equivalent to assuming 
that large-scale magnetic turbulence is homogeneous over the 
Galaxy. 



7. The results 

In this section, we present the predicted Galactic annihilation 
emission of 26 Al, 44 Ti, and 56 Ni positrons, for the various trans- 
port configurations considered. We compare these to one of 
the latest models of the allsky 511 ke V emission obtained from 
INTEGRAL/SPI observations. Last, we discuss our results in 
light of other recent theoretical studies. 

7.1. Predicted intensity distributions 

We present in Figs. [2] [3] and [4] the 511keV intensity distri- 
butions obtained in each transport configuration for 26 Al, 44 Ti, 
and 56 Ni positrons, respectively. To facilitate the comparison, 
each skymap shows the parapositronium annihilation only, for 
a positronium fraction of 0.95 and a total Galactic positron in- 
jection rate of 10 43 e + /s, of the same order as that inferred from 
observations. To emphasise the trends and differences, we show 
in Fig. [6] the longitude profiles for each skymap. To help inter- 
preting these results, we list in Table [TJ the annihilation fraction 
and the annihilation luminosities for the bulge, disk, and entire 
Galaxy. 

For a same transport scenario, the skymaps look pretty simi- 
lar: the intensity is highest in the inner Galaxy and progressively 
fades away with increasing longitude until it has dropped by 
about an order of magnitude in the outer Galaxy. The impact 
of the transport scenario on the resulting intensity distributions 
remains limited despite varying the diffusion coefficient by three 
orders of magnitude. 

In the fully collisionless scenario, positrons have very short 
ranges and they all annihilate close to their injection sites. The 
intensity distributions in that case reflect the source distribu- 
tions (and are therefore strongly driven by our assumptions about 
the latter). The normalised intensity profiles for 44 Ti and 26 Al 
positrons are almost identical, because they have a similar source 
distribution dominated by the star-forming disk. They exhibit 
a plateau over the inner +30° in longitude that actually cor- 
responds to the molecular ring at a Galactocentric radius of 
~ 4.5 kpc. In contrast, the emission from 56 Ni positrons peaks 
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at the Galactic centre since 56 Ni positrons are preferentially re- 
leased in the inner regions. It is interesting to note that when 
56 Ni positrons are confined to the vicinity of their sources (sce- 
nario A), they do not give rise to a highly -peaked 511 keV signal 
from the inner bulge, resulting for instance from massive annihi- 
lation in the strong concentration of molecular gas in the central 
~200pc; instead, they seem to annihilate in a more distributed 
fashion over the entire bulge. The results of the fully collision- 
less scenario can be expected to be representative of what is ob- 
tained when positrons are injected into the ISM with a much 
smaller average energy than that of their original /3-decay spec- 
trum, whatever the transport conditions (this is especially rele- 
vant to 56 Ni positrons that may be considerably slowed-down on 
their way out of the stellar ejecta). In that case, positrons lose 
their initially-small kinetic energy over short distances and thus 
annihilate close to their sources. 

As the diffusion efficiency is increased from scenario A to 
B and then C, the overall intensity goes down because an in- 
creasing fraction of the positrons escape the system. In the fully 
ballistic scenario, a fraction of 30-40% of the positrons manage 
to stream out of the Galaxy. The fact that the annihilation frac- 
tions are lower for 44 Ti and 56 Ni positrons than for 26 Al positrons 
is due to their higher average energy at injection, which allows 
them to survive longer and thus reach the boundaries of the sys- 
tem. This can also be seen in the longitude profiles, where the in- 
tensity drop between scenario A and B/C is stronger for 44 Ti than 
for 26 Al. The longitude profiles are not just shifted down as the 
diffusion efficiency is increased; the drop actually is stronger for 
the inner regions than for the outer Galaxy. This can be explained 
as follows: a larger diffusion coefficient allows more positrons to 
escape perpendicularly to the plane towards the halo, which low- 
ers the overall luminosity; in the same time, the most energetic 
positrons spread out radially from the innermost regions, which 
flattens the longitude profile. 

For 44 Ti and 26 Al positrons, there is almost no difference be- 
tween scenarios B and C. This is because there is almost no bulge 
contribution to the injection of these positrons, and so their emis- 
sion is mainly influenced by the transport conditions in the disk 
(which are the same in scenarios B and C). In both cases, the in- 
tensity plateau remains clearly apparent and this shows the pre- 
dominant role of the molecular ring, where positrons can slow 
down and annihilate very efficiently in the dense gas. For 56 Ni 
positrons, there is a clear change of the emission profile from 
scenario A to B and then C. In scenario B, it is interesting to 
note that despite a low diffusion coefficient in the inner regions, 
half of the bulge positrons manage to escape (see Table [T]i. A 
fraction of these are transported out to the disk and especially 
the molecular ring. As we move to scenario C, bulge positrons 
can escape much more easily towards the halo and fewer of them 
will feed the disk, the luminosity of which slightly diminishes. 
The decrease of the disk luminosity between scenarios B and C 
for 56 Ni positrons actually reflects the decrease of the molecu- 
lar ring luminosity, caused by a reduction of the supply by the 
bulge. 

The case of 56 Ni positrons reveals an important difference 
between the bulge and the disk in terms of transport. Moving 
from fully collisionless to fully ballistic transport (scenario A to 
C), the luminosity of the bulge drops by ~10% while the lumi- 
nosity of the disk decreases only by ~30% (see Table [T). It is 
obviously easier to escape the bulge than the disk. This is due to 
the relative size and emptiness of the bulge, and to the fact that 
the spatial distribution of 56 Ni positron sources in the bulge has 
a relatively high scale height of 400 pc. In contrast, the disk har- 
bours a dense molecular ring that positrons cannot easily leave, 



Table 2. Results of the fits of our predicted intensity distributions 
to about 7 years of INTEGRAL/SPI observations at 5 1 1 ke V. 



Case 


Inner 


Outer 


Disk 


MLR 


A 


1.34 ±0.17 


5.52 ±0.35 


17.1 ±1.7 (0.61) 


2687.9 


B 


1.31 ±0.17 


6.01 ±0.32 


18.8 ±1.7 (0.88) 


2687.8 


C 


1.29 ±0.17 


6.42 ±0.30 


19.3 ±1.8 (1.03) 


2684.4 


W 


1.44 ±0.17 


4.57 ±0.33 


15.6 ±1.4 


2692.3 



Note to the table: Cases A,B,C correspond to our models of all nu- 
cleosynthesis positrons with their mean injection rates and the three 
transport scenario considered. Case W correspond to the best-fit reop- 
timised Weidenspointner model. Columns 2-4 gives the fitted fluxes in 
10~ 4 ph crrr 2 s -1 for the inner bulge, outer bulge, and disk components, 
respectively. The numbers in parentheses are the scaling factors applied 
to our models in the fit. Column 6 gives the Maximum Likelihood Ratio. 



especially since most of them are injected directly in the molec- 
ular ring through massive stars, ccSNe, and prompt SNe la. The 
molecular ring is in the same time a major source and a major 
trap for nucleosynthesis positrons. 

Based on these results for each source, we computed a set 
of full emission models by adding for each transport scenario 
the contributions of all three isotopes weighted by their respec- 
tive mean positron injection rate presented in Sect. |4] To allow 
direct and meaningful comparison with the data, we rescaled 
the outputs of our transport code so that the resulting intensities 
correspond to all 51 1 keV processes (annihilation through para- 
positronium and direct annihilation of thermalised positrons) for 
a positroni um fraction of 0.97 (the value inferred from the obser- 
vations bv lJean et alJl2006h . The eventual skymaps are shown in 
Fig-El an d the corresponding longitude profiles in Fig.|7] Due to 
the preponderance of 56 Ni positrons, the intensity profile peaks 
at the position of the Galactic centre. Yet, the corresponding B/D 
luminosity remains low, with a value of order 0.2-0.4 depending 
on the transport configuration. This is an order of magnitude be- 
low the values of 2-6 inferred from the gamma-ray observations. 
In the following, we will compare our predictions to observa- 
tions and show that nucleosynthesis positrons can account for 
only a fraction of the measured annihilation signal. 

7.2. Comparison with INTEGRAL/SPI data 

From about 7 ye ars of IN TEGR AL/SPI data, w e repeated the 
analysis done in [Weidenspo intner et al. (2008b) to determine 
by model-fitting the 511 keV intensity distribution that best ac- 
counts for the gamma-ray observations. We used the same com- 
ponents adopted by these authors for their symmetric model: a 
superposition of two spheroidal Gaussian distributions for the 
bulge emission, and an exponential disk with central hole for the 
lo ngitudinally-exten ded emission (the young stellar disk model 
of [Robin e t al. 2003). All components are assumed to be centred 
at the Galactic centre. We reoptimised iteratively some param- 
eters of the model: first the widths of the Gaussians and then 
the disk scale length and height. The intensity distribution of 
our best-fit model is shown in Fig.Q] and it is si milar to the re- 
sults presented in Weidenspo intner et al. (2008b). In agreement 
with other studies, the intensity is strongly concentrated in the 
inner Galaxy with a major bulge contribution. Such a model 
seems to suggest that there is little contribution from the disk 
at \l\ > 50°, but this actually comes from an instrumental lim- 
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itation. Most INTEGRAL/SPI exposure is on the inner Galaxy, 
which is also where the surface brightness of the 511 keV signal 
is highest. Consequently, any best-fit model can be expected to 
provide a more reliable account of that part of the sky compared 
to higher longitude and/or latitude regions. Further away from 
the inner regions, the intensity predicted by the model is mainly 
driven by the assumed functional form for the fitted disk compo- 
nent. In other words, the intensity distribution at |/| > 40 - 50° 
should be considere d as presently unconstrained (see Fig. 6 in 
Bouch et et al . 2010, no significant 511 keV emission is detected 
in 5° x 5° sky pixels beyond \l\ > 40°). In the following, we will 
therefore focus on the inner |/| < 50° when comparing our pre- 
dicted intensity distributions to that derived from observations. 

Comparing this best-fit model to our predicted emission for 
all nucleosynthesis positrons, the intensity distributions appear 
to be very different. From Fig. [7] even allowing for some renor- 
malisation of the positron injection rates within the ranges of 
values presented in Sect. |4] it is obvious that the annihilation of 
nucleosynthesis positrons cannot account for the strong central 
peak observed in the inner / = +10°. So at least in the frame 
of our model, an extra source of positrons is needed to explain 
the bulge emission. Outside the bulge, however, nucleosynthe- 
sis positrons can account satisfactorily for the emission in the 
/ = +10 - 50° range, provided the positron injection rates are 
adjusted about the mean values we derived earlier (see below). 

To demonstrate this, we replaced the exponential disk by 
our predicted intensity distributions in a series of fits to the 
INTEGRAL/SPI data. The two bulge components of the best- 
fit model pictured in Fig. Q]were used to account for the bulge 
emission (only their normalisation was fitted, the morphology 
parameters were not reoptimised). The fit results are presented 
in Table[2] Our models for the disk emission can account for the 
INTEGRAL/SPI data almost as well as the exponential disk of 
the Weidenspointner best-fit model (especially if one considers 
that the bulge was not optimised for our models). The difference 
in maximum likelihood ratio between our models A or B and 
the Weidenspointner best-fit model is about 4, which would be 
an insignificant difference even if only one parameter of the disk 
model had been modified (this also confirms the above statement 
that the high-longitude emission is poorly constrained, since our 
predicted intensity profiles clearly differ from that of the expo- 
nential disk model at |/| > 50°). On the contrary, using our disk 
models seem to require a higher outer bulge emission, by ~40% 
in the most extreme case (here again, it might be necessary to 
reoptimize the outer bulge width). 

In terms of information about positron injection and trans- 
port, our fits showed that the data cannot significantly distinguish 
between the different transport scenarios tested in that work, 
most likely because our predicted intensity distributions do not 
exhibit strong variations upon modifications of the transport con- 
ditions. The best models A and B were rescaled in the fit by 0.61 
and 0.88 respectively. This could be interpreted as an indica- 
tion that the positron injection rates were slightly overestimated. 
Applying the same downscaling to all positron sources would 
imply for instance escape fractions of 61-88% and 3-4% for the 
44 Ti and 56 Ni positrons, respectively, which lie in the ranges of 
possible injection rates presented in Sect.|4] Yet, considering the 
number of uncertain parameters involved in the modelling of the 
Galactic transport and annihilation of positrons, plus the fact that 
the fit is strongly influenced by the as-yet unexplained contribu- 
tion of the bulge, such a conclusion on the positron injection rate 
should not be taken as very solid. 



7.3. Comparison with other studies 

We now compare our results with those obtained by other au- 
thors. M ost of the comparison will be done with the work of 
IHLR091 since this is one of the most extensive studies on the 
subject of nucleosynthesis positrons. 

IHLR09I claim that differential propagation of nucleosynthe- 
sis positrons can explain all properties of the INTEGRAL/SPI 
observations. We reach a different conclusion, at least for the 
morphology of the annihilation emission. We identified several 
assumptions in the work of IHLR09I that may explain the discrep- 
ancy and review them in the following. 

44 7I contribution: They assumed that 44 Ti positrons come 
mostly from SNe Ip, a peculiar subtype of SNe la. While the ex- 
act origin of the present-day s olar abundance of the daughter nu- 
cleus 44 Ca remains uncertain dThe et al.ll2006l) . we want to stress 
that the observations available today favour a dominant contri- 
bution from ccSNe (higher yield and higher occurrence rate, see 
Sect. 14.2b . hence an injection in the star-forming disk rather than 
in the central regions. 

SNe la distribution: They used a spatial distribution of SNe 
la based only on stellar mass and composed of an exponential 
disk and a bulge. We also used similar distributions (except our 
exponential disk has a central hole), but we had in addition a 
component based on star formation rate (see Sect. 14.2) . hence 
more contribution from the star-forming disk. This is somewhat 
alleviated by the fact that our stellar bulge and disk contribute 
equal rates of SNe la, while in their model the bulge is responsi- 
ble for only ~ 1/3 of all SNe la. 

Positron energies: It is not clear whether the authors used the 
complete /3-decay spectra or worked with typical energies such 
as 1 MeV or the average energies of the /J-decay spectra. This 
may have important consequences since most positrons actually 
have initial energies of a few lOOkeV only, which limits their 
range. In our study, the full yS-decay spectra were used. 

Bulge transport conditions: The authors argued that the ma- 
jority of positrons created in the inner and middle bulge within 
R < 1.5kpc are confined to that region because their scattering 
mean free path in the outer bulge is smaller and cause them to be 
reflected back. In addition, a significant fraction of the positrons 
created in the outer bulge would be preferentially transferred 
down to the inner and middle bulge. All these positrons trapped 
within R < 1.5 kpc would eventually diffuse into the cloudy 
labyrinths of the tilted disk and central molecular zone, and an- 
nihilate there. This was however not proven by a complete mod- 
elling of the entire region. 

Overall, there are many assumptions that favour a strong 
concentration of the positron creation and annihilation in the in- 
ner Galactic regions. On the other hand, our approach has its 
own restrictions, and we identified below the most important of 
them. 

The main limitation of our global diffusion model is its large- 
scale definition, with gas densities and diffusion properties being 
average values over several lOOpc at least. In that frame, it is 
not possible to simulate inhomogeneous diffusion in a realistic, 
finely-structured ISM, that is to say to follow the propagation 
of positrons through successive, identified phases with differ- 
ent transport properties. This is something Monte-Carlo simula- 
tions could do from statistical descriptions of the ISM properties 
(Alexis et al. 2012, in prep.). Regarding simply the diffusion, 
it may be that the details of the transport through various ISM 
phases of ~l-10pc sizes can be captured in large-scale diffusion 
coefficients valid over a few lOOpc scales and being possibly 
different over large Galactic regions (like the bulge and the disk 
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in our transport configuration B). Yet, in terms of energy losses, 
it might be that the average gas densities over a few lOOpc scales 
are not representative of that experienced by positrons if the dif- 
fusion really is inhomogeneous over the various ISM phases and 
positrons sample some phases more than others. This might be 
quite important especially in the inner bulge, where most anni- 
hilation seems to occur. On the other hand, the sources of nucle- 
osynthesis positrons are very likely distributed over kpc scales, 
even within the bulge. It seems improbable that diffusive trans- 
port starting from such a widespread injection could result in a 
highly-localized annihilation because of peculiar conditions in a 
few lOOpc size region. 

Another disregarded effect if the role of the large-scale mag- 
netic field topology. Whatever the nature of the transport, bal- 
listic or collisionless, the propagation of positrons is thought to 
be anisotropic on small scales, being more efficient along field 
lines than across them. We argued that the strong random fluctu- 
ations of the Galactic magnetic field over ~100pc scales would 
isotropize the propagation over large scales (the same argument 
is invoked for cosmic -ray propagation studies). Yet, the Galactic 
magnetic field is not fully random over large scales and has a 
global structure made of a toroidal part in the plane (likely fol- 
lowing a spiral pattern) and a poloidal part out of it (possibly 
X-shaped). How this global structure impact the transport of 
positrons remains unc ertain and many scenarios can be imag- 
ined (see for instance iPrantzosi 1200 6). We note, however, that 
most positrons do not travel very far from their injection sites. 
Even when they are allowed to do so (large diffusion coefficient, 
transport configuration C), about a half of the positrons annihi- 
late in sites that closely follow the initial distribution of sources 
(see TableQ}. 

Beyond the limitations of any particular model, there are 
several reasons why large-scale propagation of nucleosynthe- 
sis positrons may not occur even if interstellar conditions and 
magnetic field topology seem favourable. Most positrons actu- 
ally have initial energies of a few 100 keV only, and not 1 MeV 
as usually assumed in rough estimates, and the energy loss rate 
through ionisation/excitation and Coulomb interactions increase 
with decreasing energy below 1 MeV. This would by itself re- 
strict the range of positrons whatever their injection site. But in 
addition, all source distributions have a star-forming disk com- 
ponent, which means that many positrons are injected in the 
high-density molecular ring and annihilate there. From these ar- 
guments and the results of our modelling, it seems unlikely that 
only nucleosynthesis positrons be responsible for the observed 
annihilation emission. A corollary of the above statements is 
that the extra source needed to account for the bulge annihilation 
emission should very likely be concentrated in the inner regions, 
if the positrons have initial energies in the 100keV-l MeV range. 

8. Conclusion 

The Galaxy hosts a low-energy positron population that is re- 
vealed through a clear signature of 5 1 1 ke V and continuum anni- 
hilation coming predominantly from the inner regions. A likely 
source of such particles is the decay of radioactive species pro- 
duced by the nucleosynthetic activity of our Galaxy. We assessed 
the contribution of 26 Al, 56 Ni and 44 Ti positrons to the observed 
gamma-ray signal by simulating their transport and annihilation 
using a GALPROP-based model of high-energy particle diffu- 
sion in our Galaxy. Starting from source spatial profiles based 
on typical distributions of massive stars and supemovae, we ex- 
plored how the annihilation intensity distributions vary upon dif- 
ferent transport prescriptions ranging from low-diffusion, colli- 



sionless propagation to high-diffusion, ballistic propagation. 

In the frame of our model, these extreme transport scenarios 
at the scale of the Galaxy result in very similar intensity distri- 
butions with small underlying bulge-to-disk luminosity ratios of 
~0.2-0.4, which is an order of magnitude below the values in- 
ferred from observations. The intensity profiles are determined 
to a relatively high degree by the adopted source distributions. 
This indicates that most positrons do not travel very far from 
their injection sites, even when they are allowed to do so. We es- 
timate that at least about half of the positrons annihilate in sites 
close to their sources. As propagation efficiency is enhanced 
from fully collisionless to fully ballistic scenarios, the fraction 
of positron escaping the system increases from to 30-40%. 
The corresponding luminosity drop is more pronounced for the 
bulge than for the disk, owing to a relatively tenuous medium 
over most of the bulge volume and a high scale height of the 
bulge sources. In contrast, the disk harbours a dense molecu- 
lar ring that positrons cannot easily leave, especially since most 
of them are injected there through massive stars, ccSNe, and 
prompt SNe la. Allowing for specific transport conditions in the 
bulge does not result in a higher bulge-to-disk ratio for the anni- 
hilation emission. 

Comparing to the INTEGRAL/SPI observations of the 
511 keV emission, we conclude that the annihilation of nucle- 
osynthesis positrons cannot account for the strong central peak 
observed in the inner / = +10°. An extra source of positrons is 
needed to explain the bulge emission, and it is very likely con- 
centrated to the innermost regions if positrons have initial en- 
ergies in the lOOkeV-lMeV range. This contrasts with recent 
claims that nucleosynthesis positrons can explain all properties 
of the INTEGRAL/SPI measurements. Outside the bulge, how- 
ever, nucleosynthesis positrons can account satisfactorily for the 
extended disk-like emission, as confirmed by a fit of our pre- 
dicted intensity distributions to the data. 
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Fig. 2. Predicted 511 ke V intensity distributions for the annihi- 
lation of 26 Al positrons in each transport configuration. From 
top to bottom, the transport was assumed to be collisionless ev- 
erywhere (case A, small diffusion coefficient), collisionless in 
the bulge and ballistic out of it (case B, inhomogeneous case), 
and ballistic everywhere (case C, large diffusion coefficient). 
The given intensities correspond to parapositronium annihila- 
tion only and were normalised to a positron injection rate of 10 43 
e + /s, with a positronium fraction of 0.95. 
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Fig. 3. Predicted 511 keV intensity distributions for the annihi- Fig. 4. Predicted 511 ke V intensity distributions for the annihi- 
lation of 44 Ti positrons in each transport configuration. From lation of 56 Ni positrons in each transport configuration. From 
top to bottom, the transport was assumed to be collisionless ev- top to bottom, the transport was assumed to be collisionless ev- 
erywhere (case A, small diffusion coefficient), collisionless in erywhere (case A, small diffusion coefficient), collisionless in 
the bulge and ballistic out of it (case B, inhomogeneous case), the bulge and ballistic out of it (case B, inhomogeneous case), 
and ballistic everywhere (case C, large diffusion coefficient), and ballistic everywhere (case C, large diffusion coefficient). 
The given intensities correspond to parapositronium annihila- The given intensities correspond to parapositronium annihila- 
tion only and were normalised to a positron injection rate of 10 43 tion only and were normalised to a positron injection rate of 10 43 
e + /s, with a positronium fraction of 0.95. e + /s, with a positronium fraction of 0.95. 
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Fig. 5. Predicted 511 ke V intensity distributions for the annihila- 
tion of all nucleosynthesis positrons, using the estimated mean 
positron injection rates. From top to bottom, the transport was 
assumed to be collisionless everywhere (case A, small diffu- 
sion coefficient), collisionless in the bulge and ballistic out of 
it (case B, inhomogeneous case), and ballistic everywhere (case 
C, large diffusion coefficient). The given intensities correspond 
to the 511 keV emission from parapositronium and direct anni- 
hilation, assuming a positronium fraction of 0.97. 



Fig. 6. Longitude profiles of the 511 keV emission from 26 Al, 
^Ti, and 56 Ni positrons (top, middle, and bottom panel respec- 
tively). Black, blue, and green curves correspond to transport 
configuration A, B, and C respectively. The red curves are the 
profile of a 5 1 1 ke V sky model obtained from model-fitting to 
the INTEGRAL/SPI observations (the dashed part represent- 
ing the longitude range over which the intensity is poorly con- 
strained). The given intensities correspond to parapositronium 
annihilation only and were normalised to a positron injection 
rate of 10 43 e + /s, with a positronium fraction of 0.95. The pro- 
files were integrated over a 7r/9rad latitude band centred on the 
Galactic plane. 
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Fig. 7. Longitude profiles of the 511 keV emission of all nucle- 
osynthesis positrons for each transport configuration. The color 
coding is the same as in Fig. [6] The intensities correspond to the 
511 ke V emission from parapositronium and direct annihilation, 
for a positronium fraction of 0.97. 
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